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In the QCD context, the non-singlet structure functions of u and d -quarks are identical, save 
the initial quark densities. Electro-weak radiative corrections, being flavor-dependent, bring fur- 
ther difference between the non-singlets. This difference is calculated in the double- logarithmic 
approximation and estimated numerically. 

PACS numbers: 



I. INTRODUCTION 

Double-logarithmic (DL) contributions were discovered in Ref. [lj in the QED context and since that have become a 
popular object of theoretical investigations. On one hand, DL terms are among the most sizable radiative corrections 
in each order of the field theories at high energies. On the other hand, the ways to select the Feynman graphs 
yielding DL terms, the means to calculate DL contributions and the methods of all-order summations first developed 
in Ref. @ converted earlier examples of DL calculations into the regular technique that allows to account for DL 
radiative corrections in a quite efficient and simple way. With certain technical modifications, especially non-trivial 
for inelastic processes, the general prescriptions of calculating DL asymptotics elaborated in Ref. [2[ were generalized 
to QCD and the Standard Model of the electro- weak interactions at TeV energies where the total energy yfs >• Myy,z- 
As for the electro-weak (EW) double-logarithms, quite often in the literature they are accounted in fixed orders in 
the EW couplings. Ref. Q proved the exponentiation of the soft EW DL contributions. Such an exponentiation 
takes place for electro-weak reactions in the hard kinematics. The more involved Regge kinematics was studied in 
Refs. [3, 0]. One of the most essential difference between EW and other DL calculations is the fact that the gauge 
symmetry of the EW interactions is partly broken and the set of the EW bosons includes the massless (photons) and 
massive (W,Z) particles. The DL contributions involving soft photons are infrared-divergent and are regulated with 
the infrared cut-off /i exactly as in QED. The value of fi is fixed in final formulas with physical considerations. DL 
contributions involving soft W, Z-bosons are infrared-stable and contain, instead of fi, masses M Wl Mz of the involved 
bosons. The difference between M\y and Mz can be neglected with the DL accuracy. It makes possible to use the 
second cut-off, M (with M ^ Mw ~ Mz) instead of M w , Mz in the DL contributions from virtual W and Z -bosons. 
This approximation considerably simplifies all-order summations of EW double-logs. Another interesting topic is the 
interplay between the QCD and EW double-logarithmic contributions. In particular, it was stated in Ref. @ that the 
impact of the first-loop EW double-logarithmic terms on hadronic reactions (= EW impact) can be as large as 10% 
at energies yfs ~ 500 GeV; the role of sub-leading contributions was discussed in Ref. [2]. 

In the present paper we examine the EW impact first for the 2 — > 2 quark scattering amplitudes and then for the 
inclusive cross sections. We show that the EW impact on the elastic QCD scattering amplitudes calculated in the 
first loop is smaller than the estimates made in Ref. Q: the impact is ~ 3.5% at yfs = 1 TeV. The total resummation 
of DL contributions to the elastic scattering amplitudes increases the EW impact compared to the first-loop impact: 
the impact is 10 % at yfs = 1 TeV and growing fast with yfs reaches 30% at yfs = 10 TeV. The EW impact on the 
amplitudes of the inelastic 2 — > 2 + n -scattering of quarks can be estimated similarly. The explicit expressions for 
such amplitudes in QCD were obtained in Ref. |8| and the generalization to the electroweak processes can be found in 
Ref. [9(. On the other hand, the EW radiative corrections depend on the flavors of the involved quarks, so accounting 
for these EW corrections in DLA together with the QCD background can bring qualitatively new phenomena. For 
example, let us consider the flavor non-singlet contributions to the structure functions F± and g\ of Deep-Inelastic 
Scattering (DIS), i.e. the flavor-depended contributions to the inclusive cross sections of the DIS. They are often 
addressed as the non-singlet structure functions f^(x, Q 2 ). As is well-known, the expressions for f {±) (x,Q 2 ) include 
the coefficient functions (to describe the x -evolution), anomalous dimensions (to describe the Q 2 - evolution) and the 
initial quark densities 8q, with Sq = Su, 5d. When calculated in the QCD framework, f^(x, Q 2 ) for the u- quark and 
d- quark coincide, save difference between e^Su and e^Sd: the quark-gluon interactions do not depend of flavors of the 



2 



quarks. Electroweak corrections to bring more difference: they cause difference in the x and Q 2 -evolutions of 
the initial quarks and split into f u ^ and (the subscripts u, d label the initial quark flavors). The difference 

in the evolutions of u and d -quarks means that fu ^ even if e u Su = e 2 Sd. Impact of the electromagnetic 
~ O(o) corrections was studied in Ref. [l(| where DGLAP evolution equation[ll[ was used for accounting for the 
QCD corrections. However, DGLAP does not include resummation of the DL terms ~ a k s ln 2fc (l/a;) and the single- 
logarithmic (SL) terms ~ a* ln fc (l/a;). The point is that DGLAP was originally suggested for operating within the 
region of large x where the DL and SL contributions could easily be neglected in higher loops. Accounting for them 
to all orders in a s becomes necessary in the small- x region. DGLAP lacks the resummation, so the extrapolation of 
DGLAP into the small- a; region involves introducing the singular fits for Sq with many phcnomcnological parameters 
(see e.g. Ref. [12J) but suggests no theoretical explanations why Su and 6d should be singular. In fact, the only role 
of the singular terms in the fits is to mimic the total resummation of the leading logarithms of x(see Ref. [13| for 
more detail). When the resummation is taken into account, the singular factors should be dropped and therefore 
the fits can be simplified. On the other hand, the total resummation of the EW DL contributions to makes 
possible to estimate their impact on the small- a; behavior of the non-singlets. In doing so, we follow the approach of 
Refs. 0, H, [I3|- Through the paper we neglect the running effects for the EW couplings. 

The present paper is organized as follows: In Sect. II we consider the EW impact on the simplest and at the same 
time basic exclusive QCD process: the quark-antiquark annihilation in the hard kinematics. After that we study the 
EW impact on the inclusive cross sections in QCD. One can easily anticipate that the EW impact on the singlet 
components of the DIS structure functions should be very small: the main contributions to the singlets comes from 
the gluon ladder graphs and gluons do not participate in the EW interactions. So, we consider the EW impact 
on the non-singlet structure functions where the main contributions come from the quark ladders. In Sect. Ill we 
remind the results of Ref. 14 1 for the non-singlet structure function f^' in QCD. This expression is the solution 
of the Infrared Evolution Equation (IREE). From pedagogical reasons, in Sect. IV we first extend the QCD results 
for including the electromagnetic DL corrections. The system of IREE were all EW DL corrections are taken 

into account is obtained in Sect. V. In contrast to QCD, the evolution equations for / Uj( j involve four anomalous 
dimensions. The IREE for them are composed in Sect. VI. Besides, in order to solve the IREE for f u ,d, auxiliary 
amplitudes are obtained in Sect. VII. It makes possible to obtain explicit expressions for f u ^ first in the Mellin space 
in Sect. VIII and then in the conventional form in Sect. IX. In Sect. X we consider the small- a; asymptotics of the 
non-singlet structure functions and estimate the impact of the EW corrections on the non-singlet intercepts. Finally, 
Sect. XI is for concluding remarks. 



II. DL ELECTROWEAK CORRECTIONS TO AMPLITUDES OF 2 -> 2 SCATTERING IN THE HARD 

KINEMATICS 

In order to estimate the impact of EW DL contributions on exclusive processes in QCD, we consider the 2^2 
-scattering in the hard kinematics, i.e. in the kinematics where all Mandelstamm variables are of the same order: 

S « -t W ~U. (1) 

The impact of EW double logarithms on amplitudes of gluon scattering is surely less than the one for quarks, so we 
consider the scattering amplitude A of the annihilation of the quark-antiquark pair q q into another quark-antiquark 
pair q' q', assuming that the flavors of q and q' are different. As is known, double-logarithms are leading contributions 
among radiative corrections to this process. Calculating the DL radiative corrections to the scattering amplitudes of 
2^2 -processes in the hard kinematics is rather simple from the technical point of view because the most difficult, 
ladder Feynman graphs do not yield DL contributions in the kinematics (jTJ . For calculations with the DL accuracy 
in the hard kinematics it is convenient to use the Coulomb gauge where DL contributions arrive from the self-energy 
graphs only, which simplifies the calculations a lot. It is easy to check that this remarkable feature does not take 
place in the Regge kinematics that we consider in the next Sects. When both QCD and EW double-logarithms are 
accounted for, the first-loop contribution S is the sum of S r , with r = q, q', q, q': 

S = Y,S r = 2S u + 2S d (2) 

r 

where the subscripts u, d refer to the up- and down- quarks respectively. In the QCD context and when a s is fixed, 

S u = Sd = Sqcd = ^-a s C F ln 2 (s//i 2 ), (3) 
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with /i being the infrared cut-off. When both QCD and EW double logarithms are taken into account, 



s 



{a s C F + aQl) ln 2 (s//i 2 ) + — £ — (c' F + ^— tan 2 ln 2 (s/M 2 ) - aQ 2 r ln 2 (s/M 2 )l (4) 

sin 8w v 4 / 



where /U is the common infrared cut-off for gluons and photons. In Eqs. (|4l3p 6>vy is the Weinberg angle and Cf, C' f 
defined as (N 2 - l)/(2N) for the groups 577(3) and SU{2) respectively so that C F = 4/3, C' F = 3/4. We also use the 
conventional notations T, V and Q r for the isospin, hypercharge and electric charge of quarks. They obey the standard 
relation Q r = T3 r + Y/2, We have also used and will keep through the paper the approximation Mz ~ Mw = M. 
When a s is running, its argument in the hard kinematics is k\, so the QCD contribution (l/2)a s CF ln 2 (s//z 2 ) in 
Eqs. (|3I4[) should be replaced by 



where b is the Gell-Mann-Low function in LO. Accounting for the total resummation of DL contributions leads to 
exponentiation of the first-loop contribution S. It converts the Born amplitude Asom into 

A = A Born e- s . (6) 

Let us first estimate the impact of EW double-logarithms in the first-loop order. Defining this impact as 

e(1) = S r -S Q c D j (?) 

JQCD 

then putting Q r = 2/3 and estimating /i = 1 GeV, M — 90 GeV, we obtain that e^ 1 ) slowly grows with the total 
energy y/s: < 3% at */s < 10 3 GeV; then e^ 1 ) reaches 10% at */s ~ 10 5 GeV. Basically, fi is not fixed; in formulas 
for total cross sections it is usually replaced by suitable physical quantities like energy resolutions, minimal transverse 
momenta, etc. Here we use the estimate /i = 1 GeV because it is the scale typical for QCD and then it coincides 
with the choice of fj, for the non-singlet structure functions which we will consider in the next Sects, (see Ref. [l4| for 
detail) . Therefore, the first- loop estimate for the EW impact to the amplitudes of the 2^2 quark scattering is small 
for ^/i < 10 3 GeV. It immediately leads to the estimate < 7 % for the EW impacts on the exclusive cross sections. 
However, the situation looks more optimistic for the case of the total resummation of the double-logarithms: Let us 
define, similarly to Eq. §J), the EW impact ehard for the amplitude A in Eq. ©: 

\A-Aq CD \ 

(-hard = "j (o) 

AQCD 

where the QCD amplitude Aqcd is given by Eq. ([6]) with S = ^Sqcd- It is easy to obtain that etard grows with s 
much faster than e^, achieving w 10% at yfs ~ 10 3 GeV and exceeds 30% at \fs ~ 10 4 GeV. 



III. NON-SINGLET STRUCTURE FUNCTION AT SMALL x IN THE QCD FRAMEWORK 

The term "non-singlet structure functions" stands for flavor-dependent contributions to DIS structure functions. 
Usually, DIS structure functions are calculated with using the DGLAP evolution equations. As is known, DGLAP 
accounts for logarithms of Q 2 to all orders in the QCD coupling a s and at the same time lacks the total resummation of 
Double- and Single logarithms (DL and SL respectively) of x. Such contributions are important at small x. The total 
summation of them, including the running coupling effects, was performed in Refs. [14| with composing and solving the 
Infra- Red Evolution Equations (IREE). We will use this approach in the present paper in order to account for EW DL 
contributions, so we briefly remind below of the QCD results for the non-singlet structure functions. In order to make 
clear the fact that we discuss in this section only the QCD content of the non-singlet structure function, we will use 
the subscript "QCD" where it is necessary. Usually, notations (like /jvs) f° r the non-singlet structure functions bear 
the subscript " NS" (and the subscript "5"' is reserved for the singlet structure functions) but as through the paper 
we discuss the non-singlets only, we do not write the subscript "NS" . We denote f( +s> the non-singlet contribution to 
the unpolarized structure function F\ and use the notation f(~^ for the non-singlet contribution to the spin structure 
function gx- As is known, the latter coincides with the structure function fy. Technically, it is convenient to introduce 
the forward Compton amplitudes T^(s 7 Q 2 ) related with by the Optical theorem: 

f^(x,Q 2 ) = -QT( ± \s,Q 2 ) (9) 
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where we have used the standard notations: q is the momentum of the incoming virtual photon, p is the incoming 
quark momentum, Q 2 = —q 2 , x — Q 2 /2pq, s = {p + q) 2 ~ 2pq when x -C 1. The superscripts " ±" in Eq. © manifest 
that amplitudes have the signatures ±. It means that they are defined as follows: 

T^ = ^[T(s,Q 2 )±T(-s,Q 2 )} . (10) 

Using the signature amplitudes at high energies is absolutely necessary from the point of view of the phenomenological 
Regge theory and at the same time it is convenient technically (see e.g. Ref. [l4| for detail). Accounting for the 
summation of the DL contributions ~ (a s ln 2 (l/a;)) fe , (k = 1,...) makes necessary introducing an infrared cut-off 
p. For the sake of simplicity we identify it with the starting point of the Q 2_ev °l u tion, though it is not necessary. 
Therefore, both and depend on p as well. It is convenient (see Ref. [l4[ for detail) to use an integral 

transform to represent f^' and T^'. The Regge pole theory suggests that it should be the Sommerfeld- Watson 
transform. At s — > oo one can use its asymptotic form that looks quite similarly to the Mellin transform: 

r(±)= r^(M>) H ^ ± ) (W;2/) 



./-too 2tti \p 2 , 

where the signature factors 

£(±) _ [ e -»™ ± x y 2 ~ [! ± ! _ tm j]/2 . (12) 

As Eq. (fTTj) partly coincides with the standard Mellin transform, it is often addressed as the Mellin transform and we 
will do the same through this paper. Nevertheless, we will use the transform inverse to Eq. (|11|) in its proper form: 

2 r°° 

F (±) (cJj y) = _ / dpe-^T^ (p, y) (13) 

JO 

where we have introduced two new convenient variables p = ln(s/p, 2 ) and y — ln(Q 2 / p 2 ). Obviously, Eq. (fT3"|) does 
not coincide with the standard Mellin transform. 
Eqs. PIT)) read that 



(14) 



Evolving amplitudes with respect to p allows one to compose IREE for them. It was shown in Ref. [l4[ that in 
the QCD framework the forward Compton amplitudes T^' obey the following equation: 

T&=Tjg rn +M^®TW (15) 

where T^ rn is in the Born approximation, Mq are amplitudes of the forward quark-quark scattering. They 
should be calculated independently. After differentiating Eq. (|15[) with respect to p and applying the Mellin transform, 

Eq. (|15p converts into the following equation in terms of F^} d {uj, y): 

(u + d/dy)Fj$ D = [1 + v/2]H$ D (u)FgJ D . (16) 

The Born term T^' rn does not depend on p and vanishes after the differentiation. The term uj/2 in Eq. (fT6|) describes 
the single-logarithmic contribution. As our aim is studying DL contributions, we will neglect such SL contributions 
through the paper, though we will keep a s running. Hqq D (u>) in Eq. p6[) are related to amplitudes through 
the Mellin transform. They are new anomalous dimensions. They include the total resummation of DL and SL QCD 
contributions. IREE for Hqq D obtained in Ref. When the SL terms that do not contribute to a s are neglected, 
the IREE for H^ D is 

b {±) 

,,rr± _ °QCD , i H-(rfc) n2 , , 

where 

b QCD = a QCD + D^q CD , (18) 
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with 



a-QCD = AttA(lo)C f 



A(u,) 



dpe~ 



(p + r/) 2 



and 



dpe-^P $t[a s (s)Ta s (-s)] 



Performing integration over k 2 ± in Eq. (|20p . we obtain the following expression for D^(lu) 



QCD- 



D 



(±) 

QCD 



b 2 N J F v 



p + 77 



(p + I]) 2 + 7T 2 p + ?/ 



(19) 



(20) 



(21) 



In Eqs. (|19l2ip p = ln(s/p 2 ), 77 = ln(p 2 /Ag C£1 ) , and we have used the standard notations: C F = (N 2 - 1)/2N = 4/3 
and b is the first coefficient of the Gell-Mann-Low function. 

Eqs. (|15H21[) were obtained and discussed in detail in Ref. jl4j ]. so in the present paper we do not derive them. 



Instead, we show in next Sects how to extend the QCD results, Eqs. (|15H21[) . to the Standard Model of electroweak 
interactions. Nevertheless, let us briefly comment them. The term clqcd/{u) in Eqs. (|16I18[) is the Born contribution 
to the amplitudes of the forward quark-quark scattering, so that A{uj) is related to a s through the Mellin transform 
of Eq. (TIB")) . In contrast, the Born contribution is absent in Eq. (|16p because it does not depend on /x and therefore 
vanishes under differentiation over /j. The second term, D(uj) in Eq. (|18[) represents the approximative DL contribution 
of non-ladder Feynman graphs 1 when the s and u -channel gluons with small transverse momenta are factorized so 
that their propagators are attached to the external quark lines (see Ref. [l4j for detail). Such terms are absent in 
Eq. (|16[) because gluon propagators cannot be attached to the photon lines. The last term in the both Eqs. (|16ll7p 
corresponds to the case when a t -channel intermediate quark-antiquark pair factorizes amplitude T into a convolution 
of two on-shell amplitudes. When a s is kept fixed, A(u>) is replaced by a s and Dq^ d of Eq. (f2"Tj) is changed for 2 



D 



(±) 

QCD 



2NJ 



[IT I] 



(22) 



The relation D^ D = means that DL contributions of the non-ladder Feynman graphs cancel each other in ex- 
pressions for the forward scattering amplitudes with the positive signatures. It was first noticed in Ref. [16j in the 
QED context and remains true in QCD when a s is fixed. According to Eq. (|2ip . accounting for the running a s effects 

violates it. The expression |22|) for Dqq D (as well as Eq. (|20|) for -Dq^) consists of two factors (each in the brackets). 
The first factor (—Cf/2N) comes from simplifying the color structure t a tbt a tb of the involved graphs (t a> b are the 
S'?7(3)-gcnerators) whereas the second factor comes from integration over momenta of the virtual partons. The terms 
in squared brackets in Eq. (|2ip correspond to [a s (s) ± a s (—s)] and the exponential in this equation corresponds to 
integration of a s (k 2 _)/k 2 _. We stress that the definition of Dqgd in Eq. (f2"Tj) differs from the definition D in Ref. [ijj]: 
Dqcd = wL>. Solution to Eq. (JT7J) is 



H. 



(±) 

QCD 



uj 2 - R (± 



QCD 



(23) 



with 



B& n = 46&> n = [AttACf + ^( ± )]/(2tt 2 ) . 



J QCD — ^ U QCD 

In order to specify the general solution of Eq. (fT6j). we use (see Ref. [Uj]) the matching 

F QCD(u,y)\ v =o = F^ d (uj) , 



(24) 



(25) 



1 Through this paper we use the Feynman gauge. 

2 The sign of Eq. (31) in Ref. |14| is wrong, however this misprint does not affect the results of the paper. 
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with Fq£ D corresponding to the DIS off a nearly on-shell photon (with Q 2 = /i 2 ). It obeys the new IREE (cf Eq. (fTHJ) ) : 

"Fqcd = e 2 M") + HqcdFqcd (26) 

where e q is the electric charge of the initial quark and 5q{oj) is the initial quark density in the to -space. In contrast 
to Eq. (|16p . there is the Born contribution in the rhs of Eq. (|2T>1) because in this case we keep Q 2 ~ /i 2 , so the Born 
term depends on \i and does not vanish when differentiated with respect to \i. 

Eventually we arrive at the final answer for the non-singlet structure functions fqQ D in QCD: 

,(±> _4 r°° w ,,,,hp 



Although Eq. (|27|) is obtained for Q 2 ^> [i 2 , the shift Q 2 — > Q 2 + generalizes Eq. ([27]) to the small-Q 2 region (see 
Ref. [HI for detail). The small- a; asymptotics of fqc D is 

f$D ~ (1/^) A ^ D (28) 

where ^qcd are can ed the intercepts. Straightforwardly they can be found with applying the saddle-point method 
to Eq. (|27|) . The shorter way is to solve the equation 

^ " B$d = (29) 

for the leading singularity position and to choose its largest root. The root corresponds to the rightmost singularity 
of Eq. 427]). Ref. [3| reads that = 0.39 and A^ D = 0.42. 



IV. ELECTROMAGNETIC DL CORRECTIONS TO THE NON-SINGLET STRUCTURE FUNCTIONS 

As exchanges of virtual gluons cannot be isolated from the virtual photon exchanges, it is necessary to add the 
electromagnetic (EM) DL contributions to the QCD expression of Eq. (|27|) for the non-singlet structure functions. 
Generalization of Eq. (Tl6|) for amplitudes T^' to account for exchanges of virtual gluons and photons can be done 
in a very simple way: with replacing Hq^, d by new non-singlet anomalous dimensions h^^j accounting for both EM 
and QCD DL contributions. The IREE for h^L is similar to Eq. (fl7| : 



4mH = + {h ( EM^)f ■ (30) 



It changes Eq. (f2"7) for a quite similar expression 

A±) z\ r« 



e 2 r°° tko ,„ , , u lo „ „„(±) 



/S = 2 / ^(i/^r— brr (31) 



where new anomalous dimension sums the both QCD and EM double logarithms. It also looks like Hq^, d : 



LO - ^ 2 - 4| 

EM - 2 ' ^ ^ 

Similarly to Eq. (pM)) . is expressed through b^ M : 

4m=4m/M ■ (33) 

Now let us specify 

,(±) _,(±) +£) (±) (M) 

U EM — °QCD + «7 + ^ EM V J ^7 
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where a 7 is the electric charge of the quark: 



and 



with 



D 



(±) 

EM 



97 



47raQ 



19 



D (±) 

77 



£>(±) 
91 



4aQ 2 q C F 
b 

4aQ 2 q C F 



[1 T l]e uv 



dte-^ hit 



p{p + v) 



D (±) 

77 



4a 2 Q4[i T i] 



(35) 



(36) 



(37) 



{p + rif 



P + V 



When a s is fixed, the expressions for Dgj and become simpler 



£>(±) = £)(±) = 

97 73 



[1T1] 



(38) 



Let us explain how D^' of Eq. (|37| can be obtained from the QCD expressions for Dq^, d in Eq. (|2QII22p . Eq. IpOl 



reads that D^q^ d contains the QCD couplings depending on different arguments. 

(a) : There is cx(kjj that comes when the soft virtual gluon with momentum k 2 rs — k\ is coupled to quarks. 

(b) : There is the sum [a s (s) =p a s (— s)] from the hard virtual gluon coupled to the quarks. In and Dgy one 
of the gluons is replaced by the photon with the same momentum. In contrast to a s , we treat a as fixed: a — 1/137 . 

Therefore, when the soft gluon is replaced by the soft photon, a{k\) in Eq. (|20| should be replaced by aQ 2 and 

we arrive at . Instead, when the hard gluon is replaced, [a s (s) =p a s (— s)] should be replaced by =F 1], the 

remaining integration over k\ can easily be done and we obtain . At last, combining both previous cases leads us to 

Dffi where the both gluons are replaced by photons. This case is similar to Eq. (|22|) . save the color factor — Cp / (2N). 
Obviously, the replacements the gluons by photons change the two-gluon color factor t a tbt a tb = —Cf/(2N) for either 
t a t a = C F (for D ( ^ ] and D^) or 1 (for D$). 

In the QCD framework, the only difference between the small- a; behavior of fu^ (for up-quarks) and ffj^ (for 
down-quarks) is the difference between the initial quark densities 5u and 5d whereas both the x and Q 2 -evolutions of 
the initial up- (u) and down- (d) quark are identical, so the subscripts u and d at f^J are often dropped. Accounting 
for EM contributions brings a difference of the both evolutions on the flavor. To mark this difference, we introduce the 
non-singlet structure functions, fi and fi , with the subscripts showing the flavor of the initial quark. Obviously, 
fu^ 7^ fd even if 5u — Sd. As could be well-expected, Eq. (f3Tj) shows that the impact of EM correction on the 
small-x behavior of /( ) is very small. Indeed, the estimate of the impact cem of the EM corrections on the intercepts 
is: 



-EM 



A (+) _ A (+) 
EM ^QCD 

^QCD 



-EM 



EM ^QCD 
^QCD 



1%. 



(39) 



V. INCLUSION OF ELECTROWEAK DL CONTRIBUTIONS 



In order to include into consideration all electroweak DL contributions, adding to the gluon and photon exchanges, 
the W and Z -exchanges, we should modify the method that we used in the previous Sects, by the following reasons: 

(i) As the gauge group of the electroweak interactions is broken and electroweak bosons become massless photons 
and massive W, Z -bosons, the non-singlet structure functions acquire dependence on the both /i and M\y,z- 

(ii) MK-exchanges cause mixing of u and d -quarks, so IREE for /.„ and /lj together with IREE for the anomalous 
dimensions, are not separable (as in QCD). 

Before composing the IREE, let us introduce necessary notations. We use the notation gw for the W-coupling to 
quarks. It does not depend on the quark flavor. On the contrary, both the photon coupling e q and the Z -boson 
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coupling g q z to quarks are flavor-dependent. All these coupling are expressed through the SU(3) Standard Model 
coupling g and the Weinberg angle 9: 

9u\v = 9dw = 9w = 9/V2 , e q = g sin 9 W Q q = g sin8 w (T 3 + Y/2) , (40) 
g qZ = (g/ cos 6w)(T 3 - Q q sm 2 9 w ) = (g/ cos 9 w )(cos 2 9 W T 3 - sin 2 9 W (Y/2)) . 

We keep through the paper the standard notations T 3 ,Y and Q for the isospin, hypercharge and electric charge of 
quarks together with the standard relation Q = T3 + Y/2. We simplify the Mw,z -dependence of the non-singlets, 
assuming that in the logarithmic expressions 

M w &M Z = M . (41) 

Again, it is convenient to introduce the Compton amplitudes related to the non-singlet structure functions 

by Eq. ([9]). We will address them as the forward Compton amplitudes, although at energies y/s ^> Mw,z and 
Q 2 > Myy z the lepton and hadron participating in the DIS can exchange with 7, Z (neutral lepton currents) and W 
(charged lepton currents). In order to avoid overloading the paper we consider only the case of small Q 2 : 

Q 2 « M^ z (42) 

where the photon exchange between the lepton and quarks prevails. The other cases can be considered quite similarly. 
Under the approximation of Eq. (|41[) . the non-singlet functions f^J and the Compton amplitudes depend on 
s,Q 2 and the mass scales /i and M. We assume the following relations between the parameters s, Q 2 , M 2 , fi 2 : 

s > M 2 > Q 2 > f . (43) 
It is convenient to introduce the amplitudes F^)(u), y, z) related to amplitude similarly to Eq. 



where new variable z is introduced: z = ln(Af 2 //^ 2 ) . In accounting for DL contributions, /i acts as an infrared cut-off 
for DL terms involving soft gluons and photons whereas M acts as the second cut-off when DL terms involving soft 
W, Z -bosons are considered. In contrast to the considered above QCD and EM cases, IREE for ^ (a;, y, z) involve 

the matrix of new anomalous dimensions hit \ with i, k being = u, d, and involve the derivatives with respect to y 
and z : 

(oj + d/dy + d/dz)F& = h$(u>, z)F^ + h^ioj, z)F^\ (45) 
(w + d/dy + d/dz)F^ = hfj (u, z)F& + h<£> («, z)F^ . 

The anomalous dimensions h\t should be calculate independently. After they have been found, it is possible to find 
general solutions to Eqs. (|45|) . In order to specify them, we will use the matching 

F^(^y,z)\ v =o = Fi ± J(u J ,z) (46) 

with the amplitudes F^(lo,z). They describe the forward Compton scattering, with the EW DL corrections ac- 
counted for, in the case when the external photon has the virtuality ~ /1 2 , i.e. almost on-shell. F^J should be found 

independently (cf Eq. ([^5]) ). So, before solving Eqs. (|4"5"]) we should find h^jf* and F^J . On this step we are going 
to simplify our notations. Trough the paper we keep the DL accuracy. It gives us the right to neglect terms mixing 
amplitudes with different signatures. Therefore, all IREE we compose are separable in the signatures (see Eqs. (|16I17[) 
and Eqs. (|45l59p ). So, in what follows we basically drop the signature superscripts "(±)" but restore them when it is 
necessary. 



VI. ELECTROWEAK ANOMALOUS DIMENSIONS h lk 



In the first place let us focus on obtaining explicit expressions for hik- We will do it with obtaining and solving 
appropriate IREE. 
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A. IREE for the anomalous dimensions hik 

In our approach, in contrast to DGLAP, the anomalous dimensions can be found with composing and solving 
appropriate IREE for them. Equations for h.^ can be obtained as a generalization of Eq. (1301) : 

(uj + d/dz)h uu = b^/(&Tr 2 ) + h 2 uu + hudhdu , (w + d/dz)h ud = b™/(8n 2 ) + h uu h ud + h ud h dd , 

(w + d/dz)h du = 6f u M /(8^ 2 ) + h du h uu + h du h dd , (w + d/dz)h dd = 6f rf M /(87r 2 ) + h 2 dd + h ud h du . (47) 

The electromagnetic terms b^f and b*jj¥ in Eq. (|4"7|) are actually defined in Eq. (|M)) : 

o u « - oqcd + a utl + , o d(i - °QCD + Odd + u dd > b ud - b du — U (48) 

where 

a™ = 4™Q 2 , af/'^TraQ 2 (49) 

and D^ 1 , D dd M can similarly be taken from Eqs. (|36l37p . replacing Q q by Q u and respectively. We remind that 
we have dropped the signature superscripts " ± " for the sake of simplicity. The fact that b d AI = simplifies the 
system in Eq. (|4T|) . It is convenient to re- write Eq. (|4"7| in terms of symmetrized combinations hs,A ancl bg A A defined 
as follows: 

h h _i_ h h h h h EM h EM _1_ h EM h EM h EM h EM (Krt 

hs = n uu + h dd , h A = n U u ~ hdd , b s = b uu + b dd , b A = b uu - b dd , (50) 
and to introduce h : 

h = -u) + h s . (51) 

In these terms Eq. (|4"7| takes the simpler form: 

g = fef a V(8tt 2 ) + \h 2 + \h\ - ^ + 2h ud h du , (52) 

^ = 6r/(87r 2 ) + /^, d J^L = Kd h, ^ = h ud h. 
oz oz Oz 

Eq. Ij52")) reads that /i„d = ^d« • 

B. General expressions for hik 

Eqs. (|47l52p for ha- are non-linear, so solving them exactly is a quite serious technical problem. We do not pursue 
this aim in the present paper. Instead, we suggest an approximative procedure based on the obvious fact that the 
QCD coupling is greater than the electroweak ones. It means that in Eqs. (|47I52[) 

6f M »6f M ,6„d ,b du . (53) 

Then, Eq. (|53")) allows to conclude that 

hs > h A , (54) 



Using this relation, we can neglect /i^ and h ud h du compared to h s in the rhs of the first of equations Eqs. (|52[) and 
write an approximation for Eqs. (|52|) : 

az »7r^ 2 2 oz 8n z 

dh ud dh du 
— 7{ — = n u dh , — » — — ««d« • 
az az 
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The first of Eqs. (|55|) is the Riccatti equation and the others are linear, so they can be easily solved. The general 
solution for hs can be written as 



hs{u>, z) = uj + A 



1 + C s e> 



l-Cse 



Kd = h du = C ud exp / dth(u,t) , 



(56) 



h A 



r uEM r z 
°A 

8tt 2 



dtexp^- [ dt'h(u;,t'fj + C A exp ( dth(w,t) 



with A = y^ 2 — 2bg /(8ir 2 ) . Cs, C a (lu) and C ud (oj) being an arbitrary functions of to. They have to be specified. 
We do it, invoking the matching 



(57) 



where Hik(co) are the auxiliary anomalous dimensions corresponding to the case of the unbroken electroweak symmetry 
so that W, Z -bosons are massless, and the cut-off \i is applied to all virtual bosons. Combining Eqs. (|56|) and (|57j) . we 
express the unknown functions Cs,A,ud in terms of ■ 



where similarly to Eqs. 



Cs — -(A - H)/ (A + H) , C A — H A , C ud — H ud 

we have denoted H = -ui + Hs and Hs = H uu + H dd , H A = H uu - H d 



(58) 



C. Anomalous dimensions at the unbroken EW gauge symmetry 

IREE for Hn~ differ from Eqs. (|47| only in inhomogeneous terms: 



- *£V(80 + (H u i^ + H^H% , U H$ = b%/(8^) + H U ^H^ + fl£>fl£> , (59) 



du 



where bij^ generalize b BM to the case of the massless EW bosons. Similarly to Eq. (|34[) they can be represented as 
the sum 



',(±) 



(±) 

ik ' 



(60) 



Term 6gp D in Eq. ([60]) is defined in Eq. ([24]) , can easily be obtained from Eq. (|35j) . adding to a BM the Z and 
-boson couplings: 



a uu = add = 4naQ u + g uZ = 4tt— „ 

sin 



q (1 + Y 2 tan 2 W ) 



47ra 



^"(irf ^du 



2 2 si" 2 



sm tfvi" 



(61) 



and Djfc are generalizations of defined in Eq. (j3"6"j) . It is convenient to represent in the following way (cf 

Eq. (J3HD): 



)(±) = nW = - 



4aC f 



(1 + F 2 tan 2 6W) 



6 sin 2 6w 

oo 

dpe^P 



[lTl]e u " / dte-^lnt + 



(62) 



(3 + y 2 tan 2 6W) p(p + v) _ (1 + F 2 tan 2 6> w ) p 



=F 



4a 2 



w 2 sin 4 6*vy 
2aC* F 

~6sin 2 6> 
4q 2 

w 2 sin 4 



[iTl] 



4 (p + T]) 2 +TT 2 ' 4 

(l + r 2 tan 2 6>vi/) 2 (-l + r 2 tan 2 (V) 



16 



7_i 6 sin 0^ Jo 



up 



p 



P + V 



[IT 2] 



(-1 + Y 2 tan 2 i 



'ir) 
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When a s is fixed, the expressions for D$ and D$ look more simple and instead of Eq. ([62]) ) we obtain: 



D (±) = D (±) _ _ 

^uu — ^dd " 



8aa, 



uj 2 sin 2 aw 



Cf (3 + F 2 tan 2 ^) ri 
2 „ 4 =F -LJ 



(63) 



4q 2 (1 + Y 2 tan 2 6» 



D 



(±) 



(±) 

du 



ui 2 sin 4 0^ 
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-[1T1]- 



4a 2 



(-1 + F 2 tan 2 6> 



w ) 



u> 2 sin 4 0\v 
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[1T1] 



4a 2 



■[IT 2] 



-1 + Y 2 tan 2 %) 



Let us comment on Eqs. (|62I63[) ) . The terms ~ 1/6 in Eq. (f62|) (the term ~ a s in Eq. (|63|) ) come from interference of 
the QCD and EW DL contributions. The next term in the both Eqs. accumulate the DL contributions of the neutral 
EW bosons, 7 and Z. A part of those terms in Eq. (|62|) (all of them in Eq. (|63|) ) is proportional to the signature 
factor [1 =p 1] and therefore vanish when the signature is positive. In other words, non-ladder DL contributions to 
the amplitudes with the positive signature cancel each other totally when couplings are fixed (and cancel only partly 
when some of the couplings ar running) 3 . The presence of the last term in Eqs. (|62l63p ) demonstrates explicitly that 
accounting for the Vy-boson exchanges breaks such a compensation even when the couplings are fixed. Nevertheless, 
at fixed a s summation over flavors in Eq. (|63|) leads to the zero contribution of the non-ladder graphs: 



#) + D 

1111 1 



(+) 

ud 



D 



(+) 

dd 



(64) 



Eq. ([64]) is quite similar to the QCD result for -Dq^ with fixed a s obtained first in Ref. [13] because summation 
over flavors in Eq. ([6"3")1 is equivalent to summation over colors in QCD. As b^ are now fixed, we can solve Eqs. (|59p . 
Combining Eqs. (|60I61I62|) we see that b uu = bdd , b u d — bd u and therefore Eq. reads that H uu = Hdd and 
H u d = Hd a ■ After that Eq. ([59)) can easily be solved: 



H u 



H u 



Hdd 



1 r 
2 

bud 

~E~ 



ui — E 



(65) 



where 



7 buu 1 b u d 

Ouu ~Z o ; "ud — □ 9 ) 



E = 



\ 



46,, 



^(uj 2 -Ab uu f 



166 2 



(66) 



D. Specifying general expressions for hik 

Combining Eq. ([59"| with Eq. ([55)) and substituting them into Eq. ([56]) leads to explicit expressions for hik 



Denoting 



we obtain that 



h s (ui, z) = uj + A 



(X- E) - (\ + E)e> 



(A - E) + (A + E)e> 



h u d = h du = ^ exp f dbh(u, t) , 
E Jo 



(67) 



uEM rz 

h A = -^-5- / dtexp 



dt'h(u), t'j) exp I dth(uj,t) 



X/E = tanh/3 



A 



tanh(Az/2 + /3) 



(68) 



(69) 



3 We remind that this compensation was first noticed in Ref. [16H in the QED context. 
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Substituting it into Eq. (|67p leads to explicit expressions for hs,hA, h ud : 



A 



sinh 2 p 



h EM 



tanh(Az/2 + (3) 
1 



8tt 2 2Asinh 2 (Xz/2 + /3) 



E sinh 2 (Az/2 + /3) ' 
Xz - sinh 2/3 + sinh(Az + 2/3) 



(70) 



Eq. (|70|) manifests that the breaking the SJ7(3)(S)E/(1) symmetry of the electroweak gauge group leads to the non-zero 
Ka in contrast to the expressions Eq. obtained under the assumption of the unbroken EW symmetry. 

VII. AUXILIARY AMPLITUDES F u <d 

In the present Sect, we calculate the auxiliary amplitudes F u , F u in order to use them in Eq. (|46[) . IREE for them 
are similar to Eq. (|4l))) . save two points: the first is the absence of the y -dependence because Q 2 = for F u , d and 
the second is appearing initial contributions because they depend on \i at y = 0: 



(oj + d/dz)F u = e 2 u 5u + h uu (uj, z)F u + h ud (ui,z)F d , 
(ui + d/dz)F d = e 2 d Sd + h du (uj,z)F u + h dd (uj,z)F d . 



(71) 



The factors 8u, 5d in Eq. (|7ip stand for the initial quark densities in the u> -space. As the anomalous dimensions 
have been found in the previous Sect, (see Eq. (|TD|) ), we can solve Eq. (|71[) . Our strategy is to find a general solution 
to Eq. (|7ip and after that to specify it with using the matching to the other auxiliary amplitudes 4> u ,d of the same 
process, however obtained under the assumption of the unbroken SU(2) (8 U(l) symmetry: 



F u \z=o — 



Fd\z=0 — 



>d ■ 



(72) 



A. General solution to Eq. (|71f) 



Introducing the symmetrized combinations 

Fs = F u + F d , Fa — F u — F d , 
we can rewrite Eq. (|7ip in the symmetrical form: 

dFs/dz = {e 2 Ju + e 2 d Sd) + ( - uj + -h s (oJ, z))f s + h ud (uj, z)F s + -h A (co, z)F A 



dF A /dz= {e 2 Ju - ejSd) + - uj + -h 8 (u, z))F A - h ud (oJ, z)F A + ~h A (u), z)F s 



1 



(73) 



(74) 



It is easy to write down a general solution to Eq. (|74p in terms of integrals of h^. However, the expressions for are 
rather complicated, which makes scarcely possible performing those integrations. Instead, we obtain an approximative 
solution to Eq. ([73)1. having noticed that according to Eq. ([70]) h$ ~> hA 1 h ud . It gives us the right to drop the term 
HaFa in the first of Eq. (fTl)) . After that we arrive at the following results: 



0aH+c a H / die-**"'*) + 



<t>s 



F s = 

F A = ........ , ..... 

>0 

where cs = e^Su + e 2 d Sd, ca — e^Su — e 2 d 8d and 

1 



(75) 



dth A (w,t) + Y j dth A {uJ,t) I dve- 9 ^ 



,*(ta(,a) 



*(w,2) = J dt -uj + -h s (uj,t) =-^-ln( 



'sinh(Az/2 + /3) 
sinh/3 



(76) 



Obviously, Fs = 4>s an d F A = <Pa a t z — in accordance with the matching of Eq. (|72p. Now we should find <fis,A m 
order to specify Eq. (|75"|) . 
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B. Amplitudes 4> u .d 

Amplitudes 4> u ,d describe the forward Compton scattering off u and d -quarks under the assumption of unbroken 
EW symmetry and with the photon being on-shell. Obviously, they obey the following IREE: 



U(t>u = e 2 u 8u + H uu (j) u + H ud 4> d , 
U( t>d = z 2 d 5d + H du 4> u + H dd 4> d , 



with the obvious solution: 

(t>s = 



+ <Pd 



cs 



OA 



w — H uu — H ud 
We have used in Eq. ([78)1 that H uu = H dd . 



u — H uu + H ud 



(77) 



(78) 



C. Specifying the general solutions for F$ ,a 

When 4>a and are known, the general expressions in Eq. (|75[) can be specified: 

e -uz/2 sinh p 4 sinh ( Az / 4 ) cos h( Az /4 + /3 - ip) 



Fs = c s 



(u - Huu - H ud ) sinh(Az/2 + j3) V^ 2 - A 2 sinh(Az/2 + 0) 



Fa = ca 



z / 2 sinh 



4smh(Az/4) cosh(Az/4 + - (p) 



{uj - H uu + H ud ) sinh(Az/2 + 0) V^ 2 - A 2 sinh(Az/2 + 0) 

sinh/3 f z ,. , , , s 4 



(w — Ami — #ud) Jo 

where we have used the notation 



dth A (u),t) + 



(79) 



cs 



z/2 



V^ 2 - A 2 7 

X/uj = tanh . 



2 sinh(Az/2 + /3) 
di ^(w, t) e wt/2 sinh(At/4) cosh(A^/4 + - <p)] , 



(80) 



VIII. EXPLICIT EXPRESSIONS FOR THE ELECTROWEAK AMPLITUDES F u , d 

In the previous Sects, we obtained explicit expressions for the electroweak anomalous dimensions hik and the 
auxiliary amplitudes F Uyd . Therefore, we can now find solutions to Eq. (|4"5|) for amplitudes F u ,d ■ As Eq. (|45|) is 
quite similar to Eq. (|7ip , solving it can be done in the same way. Again it is convenient to introduce the symmetrized 
notations 



F s = F u +F d , F A = F u - F d 
and express the solution in terms of them. Obviously, 

F s (w, z-y,z) = Fs(u, z ~ I/ ) e *(«.»)-*(«.(»-i')) = 
e~ uz / 2 r sinh/3 



(81) 



(82) 



sinh(Az/2 + /3) 



w — H uu — H ua 
1 



dte^' 2 sinh(Ai/2 + /?) 



F A (u, z-y,z)= F A (u,z-y) + - dth A {Lu, t)F s (oj, z - y, t) e -*(<-,*)+*(«.*-») 

2 Jz-y 



„*(u(,«)-*(u,a-i/) 



oz/2 



CA 



sinh /3 



/^^/ 2 sinh(At/2 + ^)) + f( pi 



sinh(Az/2 + f3) L W - ff„ u + J 

+ / V dth A (t) [ due^ 2 sinh( Aw/2 + /3) + / dtfu(t) / " due" u/2 sinh(Au/2 + /?)) 



dth A (t) 



where F$ and Fa are defined in Eq. (f79|) and /i^ is given by Eq. ([70]). We remind that cs — e u Su + e d Sd and 
ca = e u 5u — e 2 d Sd. 
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IX. EXPRESSIONS FOR THE NON-SINGLET STRUCTURE FUNCTIONS 



Now we can write down explicit expressions for the non-singlet structure functions including the total resummation 
of QCD and EW double-logarithmic contributions. We express the non-singlet structure function f u of u -quark and 
the non-singlet structure function fd of d -quark in terms of their symmetrized combinations /<j and f a'- 

fs = fu + fd, /a = fu - fd • (83) 
Combining Eqs. (fl"4]) and (j82|) leads us to the following expressions: 

1 r°° dio / s \ w 

/s = 2/ joo ^(i?) Fs ^' z '"h (84) 

The Mellin amplitudes Fs,a in Eq. (|84|) are given by Eq. (|82| . When the non-singlet structure functions f u , fd are 
calculated in the QCD framework, the difference between them, ^ 0, only if ca = e 2 5it — e 2 , Sd ^ 0. Including 
the EW corrections changes the situation cardinally. Indeed, Eq. (j52"|) manifests that the expression for fA includes 
the contribution proportional to ca and, in addition, the contribution proportional to cs — e^Su + e d Sd. The latter 
contribution arises because of mixing u and d -quarks through W -boson exchanges. It means that, with the EW 
corrections accounted for, f u ^ fd even if ca = 0. We remind that Eqs. (|84|) describe /„ and fd in the region (|43|) . 



X. IMPACT OF THE EW DOUBLE-LOGARITHMS ON THE NON-SINGLET INTERCEPTS 



Let us consider the small- a: asymptotics of fi^ and /i • When they are calculated in the QCD framework, they are 
identical, save the difference between e^Su and e d 5d, and given by Eq. (f2"8")) . Accounting for the EW DL contributions 

keeps the Regge form of the asymptotics but changes the QCD intercepts Aq^ for the new ones which we denote 

A^). According to Eq. ([25]) , the intercepts are the rightmost singularities of F$,a in Eq. (j84"|) The leading singularity 
is the square root branching point in Eq. (f66|) : 

(uj 2 - 46 tm /(8^ 2 )) 2 - lQ{b ud /{SK 2 )) 2 = . (85) 

The terms b uu , b u d in Eq. (|55"|) are defined in Eq. (|60[) . They depend on the signatures, so from now on we should 
once more write explicitly the signature superscripts " ± " . It is interesting to note that Eq. (|85|) corresponds to the 
unbroken SU(3) SU(2) £g) U(l) gauge symmetry and therefore can be rewritten in the following way: 



where a 



SU{2) 



2 r 

7T 



A(u)C F + a su(2) C' F + a u{1) [Y/2)' 



2ir 2 



a/sh^Ow, au{\) = a/ cos 2 9 w ; then, C' F = 3/4, N' = 2, Y = 1/3 and 



D (±)_ D (±) c 2a lu(2) C F 
U — U QCD + 



Ac 2 Y A 



^&SU(2)CfC f 



T— — 

p + r) 



dte^' lnt 



LU 2 N' 

die-"* hit 



16a; 2 

4a [/(1) C l F y 2 
46 



-C 



8ag[/( 2 )aLf(i)C l f F 2 
4w 2 



p(p + v) 



dpe'^ff 



^ V (p + ?]) 2 + 7T 2 

p(p + r /) _ p 



(p + 7?) 2 



P + V 



(86) 



(87) 



In Eq. (Jill) we have denoted C = [1 T 1] ■ 

When a s is assumed fixed, Eq. ([56]) looks more simple: 



d(±V^ 2 = , 



(88) 



with 
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Eq. ([88)) can easily be solved analytically, the solutions, u>q are 

4 +) = , = V (« + ^a 2 + 4dF>)/2 . (90) 

On the contrary, Eq. (f86|) cannot be solved analytically. Numerical solutions to Eq. (|86|) depend on ?y and their 
maximums which we call the intercepts 4 are 



A (+) = 0.373 , A ( ~ } = 0.354 , (91) 
while the QCD intercepts ^qcd obtained in Ref. [3] are 

= °- 385 . A qcd = 0-423 . (92) 

However, the QCD intercepts of Eq. ([92| include both DL and single-logarithmic (SL) contributions. When, in 
addition to DL terms, only the SL terms contributing to a s are taken into account and other SL terms are neglected, 
the QCD non-singlet intercepts ^qcd differ from '■ 

&qcd = 0.346 , A { Q( l D = 0.389 . (93) 

Therefore, the impacts ^qcd 01? the SL QCD corrections on the non-singlet intercepts are 

A (+) _a(+) a ( - } -A m 

J+) _ ^QCD QCD _ w - _ ^QCD QCD _ „ w (Q ., 

e QCD - ~7J) 11/0 i e QCD ~ ~7T) 9/o . (94J 

QCD ^QCD 

On the other hand, impacts of the DL EW corrections on the DL QCD intercepts ^qcd are °f the same Slze: 

e (+) = y^D ^ g% e (_) = ^QCD ^ _ Q% (Q5) 

QC-D QCD 

Confronting Eq. (f3"9")) to Eq. (f9"5|) manifests that the impact of all EW DL corrections on the non-singlet intercepts 
is much greater than the impact of the electromagnetic DL corrections. It also interesting that EW DL corrections 
work opposite ways: they increase Aq^ d and decrease Aq^, which makes smaller the difference between the 
asymptotics of the non-singlets Fx and gx- A qualitative explanation to that can be easily found from considering 
Eq. (|9U|) : the expression for a in Eq. (f8"9")l manifests that adding the EW terms (all they are positive) to the QCD 

term 8a s /37r increases a and therefore increases Wq compared to its QCD value -\/8a s /37r. In contrast, there is an 
interplay between the increase of a and decrease of in the expression for cJq . Indeed, the QCD term 8a 2 /9 
in the expression for S-~^ is suppressed by the negative EW contributions (the largest of them, the second term, is 
w — 40aa s ). It means that V a 2 + 4c?(~) < a and therefore < luq + \ 



XI. CONCLUSION 



We have considered the interplay between the QCD and EW radiative corrections in the double-logarithmic ap- 
proximation. We accounted for the running QCD coupling effects but neglect them for the electroweak couplings. In 
the first place we estimated the impact of EW double logarithms on the amplitudes of the exclusive QCD processes, 
considering the annihilation of a quark-antiquark pair into another quark-antiquark pair of different flavor where the 
DL are known to be the leading contributions. We found that the EW impact for the amplitude of this process, 
though grows fast with the total energies -^/i, is less than 10% for ^/s < 10 3 GeV. As could be anticipated, the EW 
impact is more essential when the total resummation of the DL contributions has been done than for accounting for 
the double logarithms in fixed orders in the couplings. 



4 See Ref. 0, Q]| for detail. 
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On the other hand, accounting for the EW DL corrections can lead to qualitatively new phenomena which are 
absent in the QCD context. As an example, we have considered the EW impact on the non-singlet structure functions 
f^> at small x where accounting for DL contributions is known to be absolutely necessary. In order to calculate f^' 
taking into account both QCD and EW corrections in the DLA, we applied the same method of composing Infrared 
Evolution Equations that we had used for calculating in QCD. The EW couplings to quarks are sensitive to the 
quark flavors, so the Q 2 and x -evolutions of u and d -quarks are different. Besides, exchanges with virtual VF-bosons 
mix u and d -quarks. So, accounting for the EW corrections changes the QCD evolution equation of Eq. (|27p for the 
system of more involved equations in Eq. (|45p. Instead of two non-singlet anomalous dimensions Hq^, d in Eq. (|23p . 
Eq. (j4"5)) involves eight of them: , with i,k = u,d. They obey the system of non-linear differential evolution 
equations obtained in Eq. (|47f . The approximative solutions to Eq. (|47f were obtained in Eqs. (|70)) . They were 
used to obtain the explicit expressions of Eq. (|83[) for the non-singlet structure functions f u and fd in the kinematic 
region Eq. (|43|) . Besides, the expressions for H ik in Eq. (|83|) can also be used to obtain amplitudes M ik of the forward 
annihilation of quark-antiquark pairs with flavor i into the quark-antiquark pairs with flavor k: Mijr = 8tt 2 H^ . 

In the QCD context, the only difference between the non-singlet structure functions f u and fd is reduced to the 
difference in their initial densities e 2 u bu and e 2 d 3d , whereas their coefficient functions and anomalous dimensions are 
identical. In contrast, Eqs. (|83I82[) manifest that with the EW corrections taken into account, /„ - / even if 
e 2 5u — e 2 d Sd . Eqs. (|83I82|) can also be used for estimating the x and Q 2 -dependence of the asymmetry 

A i ni\ fu{x,Q 2 )~ f d (x,Q 2 ) 

AudM)= i^mrwm (96) 

in the kinematic region Eq. (|43j) . However, the discrepancy between f u and fd does not bring much difference to the 
small- x asymptotics of f u and fd '■ they both are of the Regge type with identical intercepts. Nevertheless, Eq. (f9Tj) 
demonstrates that the EW corrections change the values of the QCD non-singlet intercepts obtained in Ref. [l4| and 
reproduced in Eq. (|92|) . It is also interesting to notice that DL contributions of non-ladder Feyman graphs produce 
opposite influence on the values of the non-singlet intercepts: In the QCD framework, the intercept ^qqd °f ^ nc 
non-singlet contribution to the structure functions Fi t % is less than the intercept of the non-singlet contribution 

to g\. Eq. (f91~j) shows that accounting for the EW corrections reverses this situation. Then, Eqs. (|9H93j) manifest 
that the impact of DL EW corrections on the non-singlet intercepts is comparable with the impact of the sub-leading, 
i.e. single-logarithmic QCD contributions and reaches ~ 11%. As the intercept is the exponent in the expressions 
~ s A for the Regge asymptotics, the 11% change of the intercept due to the EW contributions is quite substantial. 
Finally, we would like to stress that similar incorporating EW corrections into the QCD expressions for the flavor 
singlet structure functions at small x should bring really small impact because the small- x behavior of the singlets is 
mostly controlled by gluon contributions. 

XII. ACKNOWLEDGEMENTS 

We are grateful to R.K. Ellis who drew our attention to the problem of interplay between strong and electroweak 
interactions. We also grateful to D.A. Ross for useful remarks concerning the EW impact on the exclusive QCD 
processes. The work is partly supported by the Russian State Grant for Scientific School RSGSS-5788.2006.2 



[1] V.V. Sudakov. Sov. Phys. JETP 3(1956)65. 

[2] V.N. Gorshkov, V.N. Gribov, GV. Frolov, L.N. Lipatov. Yad.Fiz.6(1967)129; Yad.Fiz.6(1967)361. 

[3] V.S. Fadin, L.N. Lipatov, A.D. Martin, and M. Melles, Phys.Rev. D61 (2000) 094002 

[4] A. Barroso, B.I. Ermolaev, M. Greco, S.M. Oliveira, S.I. Troyan. Phys.Rev.D69(2004)034012. 

[5] B.I. Ermolaev, M. Greco and S.I. Troyan. Phys. Rev. D67 (2003) 014017 

[6] S. Moretti, M.R. Nolten and D.A. Ross. Phys. Rev. D74 (2006) 0973 01 flhep-ph/0503152"| ). 
[7] S. Moretti, M.R. Nolten and D.A. Ross. Nucl. Phys. B 759(2006) 59 (hep-ph/0606201 ). 
[8] B.I. Ermolaev, L.N. Lipatov. Int.J.Mod.Phys.A4 (1989)3147. 
[9] B.I. Ermolaev, S.M. Oliveira, S.I. Troyan. Phys... Rev. D 66 (2002) 114018. 
[10] A.D. Martin, R.G Roberts, W.J. Stirling, R.S. Thorne. Eur. Phys. J. C39 (2005) 155. 

[11] G. Altarelli and G. Parisi, Nucl.Phys. B126 (1977) 297; V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys.lb (1972) 438; 

L.N. Lipatov, Sov. J. Nucl. Phys. 20 (1972) 95; Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 
[12] G. Altarelli, R.D . Ball, S. Forte and G. Ridolfi, Nucl. Phys. B496 (1997) 337; Acta Phys. Polon. B29(1998)1145; A. Vogt. 

|hep-ph/0408244| 



17 



[13] B.I. Ermolaev, M. Greco, S.I. Troyan. Phys.Lett.B 622 (2005) 93; hep-ph/0511343 

[14] B.I. Ermolaev, M. Greco and S.I. Troyan. Nucl.Phys.B 594 (2001)71; ibid 571(2000)137. 

[15] B.I. Ermolaev, M. Greco and S.I. Troyan. Eur. Phys. J. C50(2007)823; Eur. Phys. J. C51(2007)859. 

[16] V.G. Gorshkov, L.N. Lipatov, M.M. Nesterov . 1969. Yad.Fiz.9 (1969) 1221. 

[17] R. Kirschner, L.N. Lipatov. Nucl. Phys. B 213 (1983) 122. 



